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We report results of extensive Dynamical Monte Carlo investigations on self-assembled Equilib- 
rium Polymers (EP) without loops in good solvent. (This is thought to provide a good model of 
giant surfactant micelles.) Using a novel algorithm we are able to describe efficiently both static 
and dynamic properties of systems in which the mean chain length (L) is effectively comparable to 
that of laboratory experiments (up to 5000 monomers, even at high polymer densities). We sample 
up to scission energies of E/kBT = f5 over nearly three orders of magnitude in monomer density 
<j>, and present a detailed crossover study ranging from swollen EP chains in the dilute regime up 
to dense molten systems. Confirming recent theoretical predictions, the mean-chain length is found 
to scale as (L) cx <f> a exp(5E) where the exponents approach = 5d = 1/(1+7) ~ 0.46 and 
a s — 1/2[1 + (7 — l)/(ud — 1)] w 0.6, S s = 1/2 in the dilute and semidilute limits respectively. The 
chain length distribution is qualitatively well described in the dilute limit by the Schulz-Zimm dis- 
tribution p(s) w s 7_1 exp( — s) where the scaling variable is s — yL/{L). The very large size of these 
simulations allows also an accurate determination of the self-avoiding walk susceptibility exponent 
7 ~ 1.165 + 0.01. As chains overlap they enter the semidilute regime where the distribution becomes 
a pure exponential p(s) = exp(— s) with the scaling variable now s — L/{L). In addition to the above 
results we measure the specific heat per monomer c v . We show that the average size of the micelles, 
as measured by the end-to-end distance and the radius of gyration, follows a crossover scaling that 
is, within numerical accuracy, identical to that of conventional monodisperse quenched polymers. 
Finite-size effects are discussed in detail. 



1. Introduction. 



Systems in which polymerization takes place under con- 
dition of chemical equilibrium between the polymers and 
their respective monomers are termed "equilibrium poly- 
mers" (EP). An important example is that of surfactant 
molecules forming long flexihle cylindrical aggregates, 
so-called giant micelles (GM)EJ, which break and recom- 
bine constantly at random points along the sequence 
(see Fig. Similar systems of EP are formed by. liq- 
uid sulfurcnJ, seleniumu and some protein filamentscl In 
the surfactant literature (e.g.El) giant micelles are often 
referred to as "living polymers" although this is poten- 
tially confusing since they are distinct from systems that 
reversibly polymerize stepwise, in the presence of fixed 
number of initiators, for which this term has previously 
been reservedBQ. As direct imaging methods clearly 
demonstrate!!! GM, which behave very much like conven- 
tional polymer chainsEI, may become very long indeed 



with contour lengths up to w Ifim. However, the con- 
stant reversible scission of the chains offers an additional 
stress relaxation mechanism in comparison with conven- 
tional '!fjMn£hcd" polymers whose identity is fixed for 
all timeS&B. 



*to whom correspondence should be addressed. 
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FIG. 1. Sketch of model: Bonds of EP chains break and 
recombine constantly with rates k s — exp( — (E + B)/fcgT) 
and k r = exp(— B/k,BT) depending on the scission energy 
E and the activation energy barrier B, both supposed to be 
independent of monomer position and density. Closed rings 
and branching of chains are not allowed within our model. 

EP are intrinsically polydisperse and their molec- 
ular weight distribution (MWD) in equilibrium is 
expectectald to follow an exponential decay with chain 
length. So far, we are not aware of any direct experimen- 
tal measurements of the MWD in such systems. Of cen- 
tral interest is the mean chain length, and for GM.Lb.eiii 
has been experimental and theoretical controversyErEjEj 
concerning its dependence on volume fraction, described 
by the growth exponent (L) cx <fi a . A scaling theory 
(summarized below) gives a ~ 0.6; although this is con- 
sistent with some .data on ionic micelles at intermediate 
or high salt levelsa, a much larger exponent a ~ 1.2 is 
suggested by experiments on lecithin-in-oil reverse mi- 
celles and some nonionic aqueous surfactantsliS. Thus 
the experimental evidence concerning the equilibrium 
growth law of GM remains controversial. 

Given the shortcomings of any approximate analytical 
treatment and the difficulties with the laboratory mea- 
surements, numerical experiments, being exact within 
the framework of the respective model and able to ac- 
count explicitly for various factors that influence exper- 
iments, should help much in understanding the thermo- 
dynamic behavior and the properties, both static and dy- 
namic, of EP. However .jap to now only a small number of 
simulational studiescjEa exist, in contrast to numerical 
experiments with conventional polymers. Indeed, while 
the connectivity of polymer chains and the resulting slow 
dynamics render computer simulations a demanding task 
in its own terms, the scission-recombination processes, 
which are constantly underway in EP, impose additional 
problems for computational algorithms, mainly in terms 
of data organization and storage. Since chains constantly 
break while other fragments unite into new chains, ob- 
jects can lose their identity or gain new ones at each step 
of the simulation. 

In earlier Monte Carlo (MC) simulations on EP0 the 



systems of polydisperse polymer chains were mapped 
onto an asymmetric Potts model, in which different spin 
values were taken to represent bonded and nonbonded 
monomers as well as vacancies in a lattice. Such mod- 
els are very efficient for studying static properties of EP 
since at each update of the lattice all sites are assigned 
new spin values subject to a Boltzmann probability; dy- 
namically however this violates the connectivity of the 
chains. Accordingly this approach faithfully reproduces 
static properties, but since the kinetics of such models 
is fictitious these cannot be extended to study dynamics 
which is one of our goals. 

In a recent work by Y. Rouault and one of us, a Dy- 
namical Monte Carlo algorithm (DMC) was proposedE], 
based pu_thc highly efficient bond fluctuation model 
(BFMjEJcZl which is known to be very accurate in repro- 
ducing both static and dynamic properties (Rouse be- 
havior) of polymer chains in melts and solutionsEa. How- 
ever, the data structure used was based on a quenched 
polymer algorithm and was therefore rather slow and 
memory consuming; a radical new approach is required. 
Below we present a method to deal with these problems 
efficiently and are thereby able to study much larger sys- 
tems. 

In several EP systems, the behavior is_atrongly af- 
fected by the presence of polymeric ringsEa. For rea- 
sons that are not entirely clear however, ring-formation 
seem to be negligible in other cases, included that of 
GM. (This does not necessarily exclude a small num- 
ber of closed loops in GM systems .as may sometimes 
be seen using directimaging methodscl.) Since the latter 
are among the most widely studied examples of EP, our 
results on the ring-free case are presented here. (This 
elaborates a brief previous discussionEj.) Correspond- 
ing results on EELsystems containing rings will be pre- 
sented elsewhereEl Under sorn£LConditions GM systems 
can also contain branch pointsEj; we forbid these in the 
present work. 

After recalling some analytical predictions in Section || 
we discuss our new approach in detail in Section ||. With 
this algorithm we are able to vary the volume fraction 
over nearly three orders of magnitude and we obtain 
equilibrated systems with average chain length up to 
about (L) w 5000 (see Tab. |). We present our com- 
putational results on static properties of EP without 
rings in Section [|. A complete crossover scaling anal- 
ysis ranging from the dilute regime of swollen EP up to 
the dense Gaussian limit is performed; for systems with 
large enough chains ((£} 3> 5) wexihtain close agreement 
with recent analytical predictionsali3'li3. A delicate issue 
concerns finite-size effects which arise if one works at too 
low a temperature for any given system size. These force 
the breakdown of the (essentially) exponential MWD as 
one enters a state where a large fraction of the monomers 
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reside in a single chain. Particular care is taken in this 
regard in Section ^| In the final Section ^ we summarize 
our findings. Extension and development of these inves- 
tigations to dynamic properties of EP_.(without rings) 
will be reported in a companion paperH. 



2. Model and Analytical Predictions. 

Before describing our computational investigation we de- 
fine the physical model, and recall some essential analyt- 
ical predictions concerning various properties of equilib- 
rium conformations of EP. 

On a coarse-grained level (see Fig. |l|) systems of EP 
are characterized by the monomer volume fraction </>, the 
energy difference E between saturated and unsaturated 
bond states, the height of the activation energy barrier 
B 1 the persistence length l p and the excluded volume 
size b of the monomer (related to the cross sectional di- 
ameter of a giant micelle). Additionally, parameters for 
the non-bonded interactions may be defined. 

As approproate to the GM case we do not allow (i) 
closed rings or (ii) branching points and suppose that 
(iii) all scission energies are independent of </> and posi- 
tion of the monomer within the chain. Hence, the rates 
for scission and recombination k s = cxp(— (E + B)/kBT) 
and k r oc exp(— _B//cbT) are taken to be constant. 
For GM simplification (iii) is expected to become ac- 
curate atJaigh salinity where electrostatic effects may be 
neglectedtl 

For EP that are long compared to the persistence 
length l p , the main departure from conventional theory 
of polymer solutionst3 is that in micellar systems, the re- 
versibility of the self-assembly process ensures that the 
MWD c(L) of the worm-like polymeric species is in ther- 
mal equilibrium. This contrasts with ordinary quenched 
polymers for which the MWD is fixed a priori, and equi- 
librium only applies to the remaining (configurational) 
degrees of freedom. For EP, only the total volume frac- 
tion of the system 



SI = c(L) [In (c(L)b d ) +E + Lfi] 



(2) 
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is a conserved quantity, rather than the entire distri- 
bution c(L). It is useful to introduce the normalized 
probability distribution p{L) = ((L) / (f>) c{L) , so that 
~Y1ilP{L) = 1- The maximum chain length L maxi given 
by the total number of monomers in the system, becomes 
of relevance in the context of finite-size effects (Sec. |). 

At the level of a Flory-Huggins mean-field approxi- 
mation (MFA) the grand potential density SI for a d- 
dimensional system of EP may be written as 



where we choose energy units so that ksT = 1. The 
first term is the entropy of mixing, the second the scis- 
sion energy E of a bond and the last term entails the 
usual Lagrange multiplier for the conserved monomer 
density Eq. (Q). Without loss of generality, we have sup- 
pressed in Eq. (Q) the part of the free energy linear in 
chain length. Minimizing with respect to the MWD, 
paying attention to the constraint of Eq. ([!]), yields the 
exponential size distribution p(s)ds = exp(— s)ds where 
the chemical potential defines a scaling variable s = /iL. 
This distribution has (L}/i = 1, i.e. the scaling vari- 
able s is given in the MFA by the reduced chain length 
x = L/(L) = s. Using again Eq. (|l|) we then obtain the 
mean-chain length 



(L) = Acj) a exp(5£) 



(3) 



with a (nonuniversal) amplitude A and the MFA- 
cxponents o.mf = o~mf = 1/2. This result is expected 
to be a good approximation near the ^-temperature and 
in the melt limit where chains become free random_walks 
uncorrelated with themselves and their neighborsE3. 

Note that the chain stiffness may in principle be in- 
corporated in the above description of flexible chains by 
adding to Eq. (^|) a free energy term for the higher prob- 
ability of trans-states. It is straightforward to work out, 
that for reasonable bending energies this renormalizcs 
the scission energy E by a constant of order unity (that 
is, of order fcgT). Here we do not pursue stiffness effects 
and instead choose a persistence length l p comparable to 
the monomer size b. 

It is relatively simple to extend the above analysis to 
dilute and semi-dilute solutions of EP. We recallEa from 
standard polymer theory that the correlation length 2 
for chains of length L in the dilute limit is given by the 
size of the chain 5 = R oc L v . When the chains (at 
given number density) become so long that they start to 
overlap at L « L* w ^-Vt^- 1 ) th e correlation length of 
the chain levels off and becomes the (chain-length inde- 
pendent) size of a 'blob' 5 = £ oc L* v . Here d = 3 is the 
dimension of space, and v rj 0.588 is the swollen chain 
(self-avoiding walk) exponent. j—. 

The mean-field approach remains validE£l so long as 
the basic 'monomer' is replaced by a blob of length S. 
Across the entire concentration range we may write the 
grand potential density as 

SI = E C ( L ) [ ln (c(£)S d ) - In ((E/b) d+e ) +E + L f i 

L 

(4) 

We have taken into account here of the free energy 
change resulting from the gain in entropy when a chain 



3 



breaks so that the two new ends can explore a volume 
E d . This gain is enhanced by the fact that the excluded 
volume repulsion on scales smaller than S is reduced by 
breaking the chain; this is accounted for by the addi- 
tional exponent 8. Note that 8 = (7 — l)/v ss 0.3 with 
7 « 1.165 (as we will confirm in Sec. |4B|). In MFA 
7 = 1,9 = 0, and Eq. (g) simplifies to Eq. (g), but this 
ignores correlations arising from excluded volume effects. 
In the dilute regime Eq. (^) can be rewritten as 

Q = < L ) [Hc(L)b d ) + (7 - 1) ln(L) + E + L/i] 

(5) 

so that the relation to the well-known partition func- 
tion of self-avoiding walks (with an effective coordination 
number z), Ql oc z l L 7_1 , is recovered. 

Note that in concentrated and semidilute systems 
there are minor corrections to Eq. (|^) from the small 
fraction of chains that are too long for their excluded 
volume interactions to be screened by the surrounding 
chains (under melt conditions, this appliesE3 for those 
chains whose length exceeds (L) 2 ). This contribution is 
exponentially small and will be neglected. Note however 
that the contribution of short chains (smaller than or 
comparable to the blob size) in the semidilute regime is 
properly included in Eq. (|J) , so long as the L-dependence 
of S is taken into account. In general, this dependence 
is contained in a scaling function 



E(s) = E(s/s*) 



(6) 



Where the scaling variable s = L/j, was already intro- 
duced above, and s* = L*fj,. This function asymptotes 
to the radius of gyration and the blob size at large and 
small s respectively. 

Minimization of Eq. M) at fixed yields the MWD 



p(L)dL = p(s)ds cx ^(s) exp(— s)ds 



(7) 



which depends via S on s* . We see that the effective 
exponent 7e(s*) = s/x = (L)/i = (s) = J sp(s)ds is not 
in general equal to unity. Hence, the scaling variable s 
is not necessarily given by the reduced chain length x 
as in MFA and as supposed in ref.12-3, but by s = 7 e x. 
In the two limits far away from the crossover line one 
can readily calculate the integrals. In the semi-dilute 
limit (s* <C 1) we obtain = 7 e = 1, i.e. as in the 

MFA s = x, but in the dilute limit (s* » 1) we get 
(L)/i = 7 e = 7, i.e. s = jx. Substituting everything we 
obtain finally the distributions 



p(x)dx 



exp(-x)dx ((L) > L*) 

r^T" 1 eM~lx)dx ((£)«£*) 



We remark that any observed breaking of the scaling 
p(L)dL — p(x)dx is an indication of cither crossover be- 
tween both density regimes ((L) « L*) or finite-size ef- 
fects ((L) « 1 or (L) w L max ). One expects for semi- 
dilute configurations close to the crossover line some 
reminiscence of the dilute behavior, which should show 
up in a slightly higher j e . 

As in the MFA-case we obtain a mean chain length 



(L) = V 



(9) 



which may be cast in the generic form Eq. (g) , but with 
the exponents ay = 5d = 1/(1 + 7) « 0.46 in the dilute 
and a s = 1/2(1 + (7- l)/{vd- 1) « 0.6, 8 d = 1/2 in the 
semi-dilute limit. Thus the concentration dependences of 
the mean molecular weight (L) in the dilute and semidi- 
lute limits differ slightly form the one predicted by simple 
mean- field theory. Here finally we define (L*,4>*) quan- 
titatively as the coordinates of the intercept of the dilute 
and semidilute asymptotes on a plot of (L) vs. cf) (see 
Fig. |4A| below). Accordingly we may write 



(10) 



4>* = Pex.p(-E/ip) 
L* = Qexp(E/K) 



with exponents <p — {a s — a>d)/(8 s — Sd) = 1 + (7 + 
1)/ {vd - 1) w 3.8 and k = {vd - l)<p w 2.93. The ampli- 
tudes P and Q are similary related to the prefactors Ad 
and A s of the mean chain length (as defined in Eq. (0) 
in each density regime (see also Eq. (12) below). 

We consider finally the specific heat capacitly of the 
system. This offers a possible experimental measure of 
the typical scission energy E. Assuming this to be purely 
enthalpic and independent of temperature, the internal 
energy density U given by the density of end monomers 
U = Ej^i c(0 = Ecj)/{L). From this we get the specific 
heat per monomer 



IdU E 2 6 tt2 , m 



(11) 



(8) 



One verifies that c v has a maximum at E = 2/5 which 
shifts slightly to higher values at lower densities (from 4 
above L* to 4.4 below). 



3. The Algorithm and Configurations. 

In EP systems bonds between monomers break and re- 
combine constantly and the chains are only transient 
objects (Fig. [[J). Therefore it is relatively inefficient to 
base the data structure on the chains (such an approach 
penalized either by sorting times or waste of memory). 
Rather, one has to base it on the monomers, or even 
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better on the two saturated or unsaturated bonds of each 
monomer. This brings the algorithm as close as possible 
to what actually happens in systems of EP and makes 
it possible via pointers between bonds to avoid all sort- 
ing procedures, time consuming nested loops and arrays, 
at the expense of only one additional list required. Our 
chosen data structure is explained for three initial chains 
in Fig. [|, one of which is actually a free monomer with 
two unsaturated bonds. 




:er(-6) = i 
:er(l) = -6 
erf-5> = -1 
erf-l> = -5 
er(-2> = 4 
:er(4) = -2 



b) Recombine ibond— 2 and jbotid 



j) Break ibond—- 1 



poimer(2) = 5 



pointer)- 1) = NIL 
poiiiter(-5) = NIL 



FIG. 2. Sketch of algorithm: (a) Each monomer has two 
(saturated or unsaturated) bonds. Chains consists of sym- 
metrically connected lists of bonds: jbond = pointer (ibond) 
ipoint — pointer (jpoint). The pointers of all end-bonds 
point to NIL. (b) Recombination of two initially unsatu- 
rated bonds ibond — 2 and jbond — 5 connects the respective 
monomers imon— 2 and jmon=5. Note that only two point- 
ers have to be changed and that the remaining chains (bold 
lines) behind both monomers are not involved, (c) Break- 
ing a saturated bond ibond requires resetting the pointers of 
the two connected bonds ibond and jbond — pointer (ibond) 
to NIL. 

Using the assumption that no branching of chains is 
allowed, the two (possible) bonds of each monomer imon 
are called ibond — imon and ibond = —imon. No spe- 
cific meaning (or direction) is attached to the sign: this 
is merely a convenience for finding the monomer from 
the bond list (imon = \ibond\). Pointers are taken to 
couple independently of sign, see Fig. ||. In the pro- 
posed algorithm the bonds are coupled by means of a 
pointer list in a completely transitive fashion (jbond = 
pointer (ibond) ipoint — pointer (jpoint)) to make re- 
combinations and scissions as fast as possible; only two 
simple vector operations are required for breaking bonds 
or recombination as shown in Fig. ||. Note that this 
would be impossible in any algorithm involving only one 
bond per monomer; the latter requires an implicit se- 
quential order of segments in the chains, forcing sort- 
ing operations of the order of the mean-chain length 
for every recombination. Unsaturated bonds at chain 
ends point to NIL. Only these bonds may recombine. 
With this algorithm, no explicit distinction between end- 



monomers, free monomers or middle monomers is re- 
quired. 

As a work-horse, harnessed to this new data organi= 
zation, we exploit the bond-fluctuation model (BFM)Ea 
for the DMC simulation of polymers. This choice was 
mainly motivated by its efficiency, and the large amount 
of existing data available on monodisperse conventional 
polymers which can serve as a reference against which to 
compare the present work. However, we emphasize that 
a different choice, such as an off-lattice MC algorithmic, 
could equally well be combined with our data structure. 

In the BFM an "effective monomer" consists of an 
elementary cube whose eight sites on the hypothetical 
cubic lattice are blocked for further occupations. We 
consider the formulation of the BFM on a dual lattice, 
i.e. the center of an effective monomer is represented 
by one site-on a simple cubic lattice as introduced by 
M. MiillerEi Excluding all 26 neighboring sites from 
occupancy renders the model equivalent to— the original 
model proposed by Carmesin and KremercJ. The vol- 
ume fraction <f> is the fraction of lattice sites blocked 
by the monomers. The monomers of a polymer chain 
are connected via bond vectors b, which are taken from 
the allowed set P(2, 0, 0), P(2, 1, 0), P(2, 1, 1), P(3, 0, 0), 
and P(3, 1,0), where P stands for all permutations and 
sign combinations of coordinates. A bond corresponds 
physically to the end-to-end distance of a group of 3 — 5 
successive monomers and can therefore fluctuate within 
some range of lengths. All length are measured in units 
of the lattice spacing a. The algorithm combines typical 
advantages of the lattice MC methods with those from 
the continuous Brownian dynamics algorithm. As de- 
fined above, it corresponds to good solvent conditions 
without hydrodynamics (but respecting entanglement 
constraintsjHJ'El 

As explained above (Sec. |J) the ends of a given EP 
are not allowed to bind together in the presented study; 
before every recombination we have therefore to check 
that both monomers do not belong to the same chain. 
Because there is no direct chain information in the data 
structure this has to be done by working up the list of 
links (which adds only four lines to the source code). 
In physical time units the simulation becomes faster for 
higher E: the number of recombinations per unit time 
goes down like exp(— E), but the chain length only up 
as exp(P/2). 

The barrier energy B > is taken into account by 
setting an attempt frequency lob for scissions and re- 
combinations. This is a convenient tool for testing dy- 
namic-behavior of the system at different lifetimes of the 
chainst3, although for static properties as studied in this 
paper, the choice of B is immaterial. Therefore in al- 
most all runs reported here we set B = 0. Those sites 
of the lattice that are not occupied by monomers are 
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considered empty (vacancies) and contribute to the free 
volume of the system. We may in principle assign an 
energy — w (w > 0) for the nonbondcd interaction be- 
tween monomers in the systemEJ, and a bending energy 
s sin(fl^) with 6ij being the angle between consecutive 
bondsEZI. In the present investigation, however, we focus 
exclusively on the process of equilibrium polymerization 
of entirely flexible chains in an athermal solvent setting 
w = s = 0. 

Time is measured, as usual, in Monte Carlo steps 
(MCS) per monomer. Each MCS is organized as follows: 

• A monomer is chosen at random and allowed Jxi 
perform a move according to the BFM algorithrrO. 

• With a frequency u>b = ex p(~ B), i.e. every 1/u>b 
MCS, one of the bonds is chosen at random (re- 
member that there are twice as many bonds as 
monomers). If one of the bonds happens to be 
a saturated P(2, 0, 0)-bond an attempt is made 
to break it, otherwise if it is unsaturated, i.e. 
the monomer is at the end of a chain or a free 
monomer, an attempt is made to create a bond 
with another monomer that might be present on 
any of the 6 neighboring £(2, 0, 0) sites. Applying 
the Metropolis algorithmES a scission is performed 
whenever the value of a random number between 
and 1 is smaller than exp(— E). On the other hand, 
with the bond energy being positive (E > 0), re- 
combination is always accepted so long as a ring is 
not thereby created. 

Note that P(2, 0, 0)-bonds are broken irrespective of 
which particular bond vector, i.e. which of the 6 possi- 
bilities, they stand for, and that, therefore, detailed bal- 
ance requires that for recombination all (and not just one 
as in refBj) of these sites have to be checked for possible 
unsaturated bonds. These neighboring sites have to be 
checked randomly; a typewriterlike search along the list 
of possible bond vectors creates correlations in violation 
of the detailed balance requirement. (This has subtle 
but measurable consequences.) Note finally that our de- 
cision to restrict the breaking and recombination to the 
6 shortest (P(2,0, 0)) bonds avoids ergodicity problems 
arising from crossecL-Bairs of bonds that can result in 
immobile monomers H3 

In the presented study the volume fraction was varied 
over nearly three orders of magnitude from <p = 0.001 
{L,nax — 1) 000 monomers per box) to <j> = 0.6 (L max — 
75, 000 monomers per box) as shown in Fig. 0. For den- 
sities smaller than <j> = 0.1 cubic lattices with linear 
dimension Sb = 200 were used, for higher densities a 
smallex_box with Sb = 100. This should be compared 
to ref.E2l where a Sb = 30-box with L max — 1, 300 par- 
ticles (4> = 0.4) was the largest achievable. One should 



bear in mind that densities around 0.5 correspond to 
extremely dense systems (melt conditions) in the BFM, 
since at higher densities the blocking of neighboring sites 
by other monomers leaves no room for mov£m£nt and the 
system goes effectively into a glassy stateE§E3. 
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FIG. 3. Simulation parameters (E, <f>) in relation to the 
crossover density (f>* . Above and to the right of the dashed 
line is the semi-dilute regime, below and to the left the di- 
lute regime. The crossover density follows an exponential 
decay with (f>* » Pexp(— E/<p) where if — 3.82 and P — 0.26 
consistent with Eqs. ( ]To| ) and (|l^). At volume fractions of 
about 4> = 0.5 the correlation length becomes of order of the 
monomer size (the melt regime) . The dotted line corresponds 
to (L) — 5; to the left of this line chains are too short for good 
scaling. For some systems (filled symbols) we have in addition 
systematically varied the frequency of scission/recombination 
uib by varying the activation energy Ref£3. For certain pa- 
rameter choices (circles) we have checked carefully for finite- 
size effects. 

The starting configuration consists of randomly dis- 
tributed and nonbondcd monomers which we cool down 
step by step (a sequence of so-called 'T- Jumps') each 
step sampling a higher scission energy up to the maxi- 
mum E — 15. As mentioned above, for the static results 
we usually set B — 0, i.e. ujb — 1. Due to the constant 
breaking and recombining of the bonds the equilibra- 
tion is then much faster than in systems of quenched 
polymers. Following Rouault et alEJ we have checked 
this explicitly by monitoring the evolution of mean-chain 
length (L(t)) and radius of gyration (i?^(<)) after every 
T-Jump. 

Measurements of static properties such as the mean- 
square end-to-end distance the radius of gyration 
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(Rg) or the specific heat c„ (see Tab. |l|) were performed 
in intervals of roughly it/ 10, where r t is the terminal 
relaxation lime of the system for the parameters (</>, E) 
considered^. Typical runs (with B = 0) covered eas- 
ily up to 10 — 100 r t . The melt density (j> — 0.5 was_a 
particular focus for our study of dynamical propertiesE3, 
requiring better statistics. Hence, for this <f> we have 
sampled over 128 independent configurations. However, 
for all static properties other than the specific heat c v 
and, for the largest E values some MWD data, a long 
run from a single starting configuration was quite suffi- 
cient. Indeed, statistical errors are generally within the 
symbol size in the data presented below. 

There are two different kinds of "finite-size effects" in 
this simulation. The first, and most important, arises 
from chains that are too small. Not surprisingly, we 
found in configurations with (L) < 5 (below the dotted 
line in fig. ^) non-universal behavior. For clarity these 
data points are omitted in most of our plots. Second, 
for systems with (essentially) exponential MWD one has 
to worry about the system size whenever (L) becomes of 
the order of the total number of monomers in the box 
L m ax- This important issue will be addressed in Sec- 
tion ^ where we present results of a systematic finite-size 
study summarized in Tab. |^ confirming unambiguously 
that the systems reported are indeed sufficiently large. 

Most of the simulations on static properties were per- 
formed on two single DEC Alpha workstation over a pe- 
riod of year. A parallel version of the algorithm, how- 
ever, was also developed and some of the computations 
at melt density mentioned above have been carried out 
on the facilities of EPCC, Edinburgh. The latter compu- 
tations focused mainly on the dynamics of equilibrium 
polymers as we report elsewhereO. 

4. Computational confirmation of scaling predic- 
tions. 



In what follows we examine the influence of density 4> 
and scission energy E on mean chain length (L), MWD 
c(L), specific heat per monomer c v and the size of the 
chains, as measured e.g. by the radius of gyration (-Rg). 



A. Mean Chain Length. 



In Fig. |A|a we show in a semi-log plot the measured vari- 
ation of the mean chain length {L} versus the bond en- 
ergy E. Note that we have been able with this new algo- 
rithm to obtain mean chain length of up to (L) ss 5, 100 
(see Tab. |l]) which is comparable with (at least some) 
experimental systems of EP. Giant micelles, for exam- 



ple, are somewhat rigid and the persistence length rela- 
tively large, l p ~ 16 nm whereas R g 100 nm, result- 
ing in around 100 statistical segmentalU, although miidp 
longer chains could arise in some semidilute systemsEa. 
As mentioned in the Sec. || and elaborated further in 
Sec. H below, it is the finite size of the systems, rather 
than the equilibration time, that prevents us studying 
higher E (lower temperatures) since the largest chain 
would then comprise too high a proportion of the total 
available monomers in the system. 



For long enough chains (above the dotted line in 
Fig. |^) the chain length increases exponentially with 
scission energy E, as predicted by Eq. (||). The data 
confirm with high precision the predicted exponents 
5 d w 0.46 ± 0.01 and 5 S = 0.5 ± 0.005 in the dilute 
and semidilute regimes respectively. The growth expo- 
nents a are most readily confirmed (following the scaling 
prediction Eq. @) by directly plotting the "number of 
blobs" I = (L)/L* against the reduced density (/)/(/>*■ 
The dat a collapse onto a single master curve, seen in 
Fig. 4Ab, is indeed remarkable and is one of the main 
results of this work. The two indicated slopes are com- 
parisons with the two asymptotically predicted growth 
exponents (ay — 0.46, a s « 0.6). (Note that in the 
dilute limit finite-size effects are visible for low E val- 
ues where chains are extremely. short ((L) < 5).) This 
finding is at variance with ref.E2l where a much stronger 
growth with density is reported for systems containing 
only 8400 monomers with mean-chain lengths up to 210. 
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FIG. 4. (a) The average chain length (L) for a wide range 
of densities <j> and energies E. (b) Rescaled average chain 
length I = {L)/L* versus reduced density cj>/cj>* confirming 
the scaling Eq. (^) and the amplitudes Eq. ([12]). The two 
slopes are comparisons with the predicted growth exponents 
ad ~ 0.46 and a s ~ 0.6. 



The distribution of chain lengths at equilibrium, c(L), 
is presented in Fig. 4Ba on semi-log axes for various E 
at high density (j> — 0.5. The fluctuations in the sam- 
pled lengths increase considerably for very long chains 
where correlations between successive configurations de- 
teriorate the statistics. To the available accuracy, the 
normalized distribution p(x) plotted versus the reduced 
chain length x = Lj (L) collapse perfectly on s ingle 'mas- 
ter' curve as shown in the insert of Fig. 4Ba; thus the 
mean chain-length {L) contains all energy information. 
The exponential decay confirms Eq. (|^a). For compar- 
ison we have indicated the exp(— jx) behavior which is 



clearly not compati 
agreement with rcf.' 2 ^ 



•\e with the data. This finding is hi 
but in clear contrast to Gujratitfl 
according to whom the Schulz distribution holds inde- 
pendently of the overlap. 



By plotting A — (L) /(f) a / exp(5E) versus (L) or E one 
fits for each density regime an amplitude Ad « 3.6 ± 0.2 
and A s w 4.4 ± 0.1 respectively. From these pref actors 
the amplitudes governing <fr* and L* defined in Eq. ( fTo| ) 
can be obtained 



P = {A d /A s ) l/{(Xs ~ ad) ra 0.3 ± 0.1 
Q = A s P a ° w 1.9 ±0.5 



(12) 



(The high error bars are caused by the small difference 
in the two growth exponents resulting in the large ex- 
ponent l/(a s — ad)-) Although these are defin ed b y 
the crossing point of the two asymptotes in Fig. 4Ab, 



it is notable that the asymptotes are followed, within 
numerical accuracy, all the way to the crossing point. 
(A large deviation is of course not possible since the 
two asymptotic slopes are not very different.) Thus the 
crossover between the two regimes occurs rather sharply 
at (L)/L* = 4>/<f>* = 1, and this can be used to define a 
crossover line on the <j>, E o r (L), E plot. Indeed dashed 
lines in fig. (||) and fig. (4A 1) indicate the position of the 
crossover lines (f>* and L* using the exponents tp = 3.84 
and k = 2.93 together with the above amplitudes P and 



B. Molecular Weight Distribution. 

We consider first the MWD in the dilute (I < 1) and 
semidilute (I ^> 1) limits, far away from the crossover. 
Thereafter we discuss the crossover effects near I « 1 
where a significant fraction of the chains are still smaller 
than the blob size of the semidilute network. 



While at high densities we observe perfect exponen- 
tial scaling of p(x), at lower dilute densities (with suffi- 
ciently long chains) our results are qualitatively consis- 
tent with the Schulz distribution Eq. ( [4 Efo ) . We com- 



pare in Fig. 4Bb this prediction (bold line) with data 
sets for configurations in the dilute limit. To stress the 
systematic difference, we have included the high density 
prediction. 



As shown in the insert of Fig. 4Bb our MWD at di- 
lute densities are also qualitatively consistent with the 
additional power-law dependence p oc x J ~ 1 in the limit 
of small x. Note that the maximum of the distribution is 
at xm = (7 — l)/7 ~ 0.1 corresponding to a chain length 
L = 0.1(L) w 3.1 for = 0.005 and « 4.2 for (j) = 0.01. 
Hence, we could not expect to reproduce accurately the 
power-law regime for x ^ xm- This would require con- 
figurations of at least (L) « 1000 in the dilute regime; 
due to the finite-size effects in the range of scission en- 
ergies used (discussed in Sec. || below) this is at present 
not feasible. (Since in 2d one has 7 = 43/32 which is 
much farer from the MEA value 7 = 1 this part of the 
distribution can be more efficiently probed in 2d, as done 
recently by Rouault and MilchevEa.) Qualitatively, how- 
ever, we believe that this result is unambiguous. Note 
that we see no evidence for a possible negative exponent 
in the power law in Eq.(|8|), as postulated in soiae treat- 
ments of the unusual diffusive behavior in GMli£l. 
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FIG. 5. Molecular Weight Distribution, (a) c(L) at high 
(melt) density <f> = 0.5 for energies E as indicated in the 
figure. Insert: Data collapse of the normalized distribution 
p(L)dL — p(x)dx = exp(— x)dx. (b) Dilute limit confirming 
p(x) oc exp(— fx). Insert: MWD for small x for two config- 
urations with <j) — 0.005 (squares) and <f> = 0.01 (diamonds) 
both at E — 10 in the dilute regime. We compare (bold line) 
with the prediction of Eq. (|8|b). 

At intermediate densities, slightly above the crossover 
line, a non-negligible fraction of chains are smaller than 
the blob size, and are thus fully swollen (these chains 
may fit among the network of chains of average size 
(L) without being seriously perturbed by the interchain 
interaction). The distribution therefore crosses over 
smoothly from the dilute limit as depicted in Fig. [4 B| b 
to the semi-dilute limit of Fig. 4Bb (data not shown). 
Throughout the parameter range we have fitted system- 
atically the effective exponent 7 e from the exp(— 7 e a;)-tail 



of the MWD. The values obtained are plotted in Fig. 4B 
versus the number of blobs I. We confirm 7 e — > 7 = 1.165 
and 7 e — > 7 = 1 in the dilute and semi-dilute limit re- 
spectively. In between we observe a crossover for the 
effective exponent. Note that the error bars are mostly 
around ±0.01, however, a much higher accuracy is in 
principle feasible by this method. Note that the value 



obtained in the dilute limit compares well with the best 
renormalization group estimate 7 = 1.1615 ± 0.0011E3. 
We believe this is the most reliable simulation determi- 
nation so far of this well-known polymer exponent. 

In passing, we recall that the polydispersity index 
/ = (L 2 ) I '(L) 2 , i.e. the ratio of weight average and 
number average molecular weights, becomes I = 1 + 1/7 
and I = 2 in the dilute and scmidilutc limits respec- 
tively. This offers (in principle) an additional method to 
check the 7-exponent. However, due to the difficulty to 
measure accurately the distribution for small x (which 
contributes strongly to I), we obtain with this method 
values that are slightly larger than that quoted above for 
the dilute limit. 

To summarize, the behavior found in a large range 
of / is in support of recent treatments of the prob- 
lem by rnpfl.ns of renormalization group and pealing 
analysesEj£3't3 and in contrast to earlier claimsa that 
the Schulz distribution, Eq.(|^), will hold independent of 
the degree of overlap between the chains. 
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FIG. 6. Effective exponent 7 e versus the number of blobs 
per chain I = {L)/L*. 

C. Specific Heat. 

A quantity that may be accessible experimentally is the 
specific heat c„(see also Tab. |l|). Due to its variable num- 
ber of broken bonds (chain ends) EP should absorb or 
release energy as the temperature, i.e. E, is varied. In 



Fig. 4 C we plot the specific heat per monomer c v <f> s ver- 



sus scission energy E for two densities within the semidi- 
lute regime. We compare with the prediction of Eq. ( pi] ) 
where we have used the exponents a s = 0.6 and 8 S = 0.5 
and the amplitude A s = 4.4 estimated in Sec. [4 A| . We 
find a qualitatively good agreement, especially, as es- 
pected, for larger chains (E > 5). 

In the insert we check the prediction c v {L)/E 2 = 6 
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(plotted versus the numbers of blobs I) for all configura- 
tions with (L) > 5. Although the statistics is too poor to 
separate the extremely small difference between Sd and 
6 S the data 'collapse' around 6 = 0.5±0.1 is qualitatively 
satisfactory for this ambitious test. 

Note that although this heat capacity has some re- 
semblence to that predicted for stepwise living polymers 
in the present of initiators, the latter systems are fun- 
damentally different in that they can show a (possibly 
rounded) phase transition at nonzero temperaturetffl. In 
the present system no critical phenomena are observed 
at finite temperature; this is confirmed by the fact that 
the measured heat capacity becomes system-size inde- 
pendent in the limit of large systems (compare Tab. @). 



and recombination events 
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FIG. 7. The specific heat per monomer c„. We plot c v (f) aB 
versus scission energy E which is here equivalent to plot 
against the inverse temperature. The line is the prediction in 
the semidilute regime. Insert: Scaling-plot of c v (L)/E 2 = S 
versus number of blobs I in agreement with Eq.dlll). Symbols 



for the various densities as in Fig. 4Aa 



D. Conformational Properties. 

Following ref.B the average over all chains of the mean- 
square end-to-end distance (-Re) an d the radius of gyra- 
tion (Rg) have been measured and are listed in Tab. @. 
To obtain meaningful results for chain size, when large 
chains are present, we must however 'unwrap' the pe- 
riodic box starting from the position of one chain end, 
using the bond vectors between consequent monomers 
of the chain. Not surprisingly the mean bond length 
(fr 2 ) is nearly identical to the ones obtained in the study 
by Paul et alBa and is likewise slightly decreasing with 
density. Also we observe (for sufficiently long chains) a 
persistence length l p ss 1.16. We include in the table the 
mean interchain distance H = (8(L) /(f)) 1 / 3 which is of 
relevance in the context of diffusion controlled scission 



FIG. 8. Crossover scaling plot of the reduced mean-square 
end-to-end distance u e = (R^/Rq 2 and the reduced ra- 
dius of gyration u g — (^R 2 ^ / Rq 2 vs the scaling variable 
v — (Ro/H) 3 . Same symbols as in Fig. 4 A a, asterisks denote 
monodisperse chains. 



In Fig. 4D we demonstrate that the mean chain sizes 
for EP follow the same universal function as conven- 
tional quenched polymers. We want to compare sizes 
of chains of given mean-chain length with the size of 
swollen dilute chains of the same length; therefore we 
define Rq = b{L) u . We plot the reduced average chain 
size u e — (i? 2 )/i?o 2 and u g = (Rg)/Ro 2 for the end- 
to-end distance and the radius of gyration respectively 
as function of the scaling variable v — (Ro/H) 3 . This 
choice of variable, rather than the alternative (/)/</)* (to 
which v is proportional) enables a comparison between 
EP and conventional monodisperse polymers (asterisks 
in Fig. 4D). Data for the latter is taken from Paul 
et alc3. In the dilute regime we have R e ~ R g ~ Rq 
and the scaling function u approaches a constant as 

Note that the plateau for u e 



can be seen in Fig. 4D. 
is slightly above 1 in agreement with the persistence 
length given above. In the semi-dilute limit the chains 
are Gaussian on length scales larger than the blob size 

£ (i.e. Rg (T\^~^'^\ imnlvino' t.Vip Qpfllincr 11 rY 



(L) 



implying the scaling u oc 



v 

with exponent /3 = (2v - l)/(3zv - 1) w 0.23. It is 
remarkable that, on this scaling plot, EP and conven- 
tional monodisperse polymers are nearly indistinguish- 
able; the two universal functions are virtually identical 
to within numerical accuracy. Note the location of the 
crossover density at u* 3 ± 1, i.e. R e w Rq ss 1AH. 
A consistency check with the estimates of the ampli- 
tudes P and Q from Eq. ([l2]) gives a slightly lower value 
v* = PQ 3l/ ~ 1 (b/a) 3 /8 w 1. However, in view of the error 
in locating the crossover values and the large exponents 
involved in the estimation of P and Q this is reasonably 
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consistent. 



This conclusion is corroborated in Fig. 4D where 
we show the distribution of chain sizes, {Kg) and 
(-R 2 ) L , averaged not over all chains present (as con- 
sidered above) but over all chains of length L, plotted 
against L (rather (L)). For the melt density <fi = 0.5 
we find (Rl) L ~ 6(Rg) L oc L 2l> with a Flory expo- 
nent v « 1/2. For the much smaller semi-dilute den- 
sity tj) = 0-01 we obtain swollen chains with an excluded 
volume exponent v w 0.6 (or perhaps slightly larger). 



Gaussian {R 2 g ) L oc L. We may measure the values of L* 
and £ directly from the crossover between both regimes 
and compare them with the values obtained above. From 



the intersection point of the two slopes in Fig. 4Db one 
has £ 2 « 420 a 2 and L* « 160. Remembering that 
every monomer occupies 8 sites of the lattice one ob- 
tains as a consistency check roughly <j> = 8L* /£ 3 rj 0.13 
whi ch n early matches the actual density of the system. 
Fig. 4Eb confirms (within numerical accuracy) the sharp 
crossover betwee n as ymptotic forms as discussed in con- 
nection with Fig. 4 Ab and is another of our main results. 



5. Finite-Size Effects. 
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FIG. 9. Variation of chain size with chain length L. (a) At 
very low density (e.g. (f> — 0.05 (dashed lines)) the chain are 
swollen R oc V with v « 0.6. At density cj> = 0.5 one has a 
melt of Gaussian chains (u — 0.5). (b) For a semi-dilute sys- 
tem, as for cj> = 0.1 at E = 10 and E = 15, one can identify 
two regimes. Up to sizes £ the chains are swollen, on sizes 
larger than £ they are Gaussian. 



Fig. 4 Do shows the distribution of (i? 2 ) at an in- 



termediate density which can be used to determine the 
average size, £, of a blob containing L* monomers. For 
L«i' and (i? 2 ) <C £ 2 , i.e. within the blob, the ex- 
cluded volume interactions dominate and (R%) L oc L 12 . 
For larger L (i.e. also larger distances) the chains become 



Since no chain can be larger than the total number of 
monomers present in the system, the exponential MWD 
must break down whenever the mean chain length be- 
comes too large, i.e. when the average number of chains 
per box (M) « L max /(L) of order one. (For the same 
reason it must always break down in the high molecu- 
lar weight tail, in any case.) Note that for the highest 
energy E — 15 used we find typically an average chain 
number (M) w 10. In fact, the (rather noisy) MWD 
show even then qualitatively an exponential decay. 

In order to understand finite-size effects quantita- 
tively we have explictly performed a systematic finite- 
size study for E = 10 and E — 15 at density tfi = 0.5 
which is summar ized in Tab. [§ In Fig. |a the mean chain 
length (L) and the end-to-end distance R e reduced by 
their asymptotic values for infinite systems (taken from 
the predictions in Section |^ together with the amplitudes 
obtained above), are plotted against an obvious finite- 
size scaling variable, which is the number of monomers 
in the box L max divided by the average chain length of an 
infinite system, (L). We see that the systems for which 
we have presented results above are indeed lying in the 
asymptotic plateau region. (Additionally, this confirms 
the amplitudes and exponents presented already above.) 
The scaling curves for the mean chain length and the spe- 
cific heat (not shown) are - not surprisingly - similar, the 
specific heat being the noisier quantity. Both decay in 
the small system limit with the system size (slope one), 
implying the formation of a single long chain of length 
~ L rnax oc <fi. Hence, the effective growth exponent a is 
likely to be overestimated in computational studies on 
relatively small systems. From L max /(L) oc we 
see that increasing the density at fixed scission energy 
should decrease the finite-size effects shifting the system 
to the right in fig. [|. Note that at intermediate sys- 
tem sizes above L max / (L) « 1 the mean chain length is 
slightly larger than the asymptotic value (see Tab. |^). 

For the chain size in the limit of strong finite size 
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effects our results are qualitatively consistent with a 
square root dependency, R oc L max 1 ^ 2 1 rather than a 
1/3 exponent for a single compact chain; this is because 
the chain size is defined by 'unwrapping' the periodic 
box starting from one chain end (as mentioned above) 
and nothing prevents a single Gaussian chain from wrap- 
ping repeatedly to fill the periodic box. In Fig. ||b we see 
directly from the MWD how the transition to the single 
chain 'phase' occurs. For small systems {Sb < 30) we 
find clear peaks in the distribution at L — L max which 
disappear as we further increase the system size. This 
finite-sizes study confirms unambiguously that the con- 
figurations presented above are indeed large enough - at 
least for the static properties. It is however not so clear 
that this is still holds for dynamic properties in the limit 
of large barrier B, and to will address this issue in the 
second part of this workE3. 
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FIG. 10. Systematic variation of the system dimension 
Sb- (a) Reduced chain length (L) and end-to-end size R e 
versus L ma x reduced by the average chain length of an infi- 
nite system. Solid symbols E — 15, open symbols, E = 10. 
(b) MWD for <f> = 0.5 at scission energy E — 10 for various 
system sizes as indicated in the figure. 



6. Discussion. 

In the present Monte Carlo simulation of static (and 
dynamicE3) bulk properties of EP we have used a greatly 
improved version of a recently proposed algorithmE3. 
The original data structure based on the BFM polymer 
chain model was completely altered, recognizing that the 
fundamental entities in EP systems are not chains but 
monomers, or (in our algorithm) the bonds connecting 
these. This makes the algorithm much more powerful, 
and allows us to study systems across three decades in 
monomer density with up to 75, 000 monomers per con- 
figuration (for cb = 0.6 on a DEC Alpha workstation). 
With this algorithm we can equilibrate systems with 
mean-chain length up to (L) w 5000 in the melt, which 



is at least one 
previous stud 



s of magnitude larger than any other 
. This was achieved with negligible 
finite size effects as we have explictly checked. 

Therefore, we are able to compare simulational results 
to data extracted from laboratory observations and to 
the scaling predictions of theoretical treatments, known 
to hold in the asymptotic limit of sufficiently long chains. 
Some earlier numerical observations which showed dis- 
crepancies with analytic results could thereby be traced 
to finite-size effects arising from small systems and/or 
short chains. 

In the present investigation we simulated EP where 
the formation of rings is not allowed. (A study for sys- 
tems containing rings is underwayE3.) The nonbonded 
monomer interactions have been set to zero (athermal 
conditions) so that the overall picture was not compli- 
cated by phase separation at low temperature. For sim- 
plicity the chains are modelled as totally flexible. ._. 

At variance to some recent computational studiesEj 
we find that the static properties in the asymptotic 
limit of large (L) agree— weII with analytical predic- 
tions nf the scaling theoryEJu and renormalization group 
studieslljll3: 



1. The mean chain length varies as (L) = 
L* (<fi/(f)*) a oc <p a e-np(SE) with density </> pajtLend- 
cap energy E. As predicted analyticallyMl!3, we 
find the exponents a<j = 5<j = 1/(1 + 7) = 0.46 
in the dilute regime and ay = 0.6, Sd = 0.5 in the 
semi-dilute limit. We confirmed the scaling behav- 
ior of the mean number of blobs I = (L) / L* versus 
the reduced density </>/</>*. This shows a rather 
sharp crossover which enables clean crossover lines 
cb*(E) and L*(E) have been located. 

2. In the dilute limit the MWD scales consistently 
with a Schulz distribution p(s) oc s 7 ~ 1 exp(— s) 
with a scaling variable s — jL/(L). In the 
semi-dilute regime for large enough (L) the MWD 
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decays exponentially with chain-length p(s) = 
exp(— s) where the scaling variable becomes the 
reduced chain length s — L/(L). Between both 
limiting regimes we observe a relatively gradual 
crossover at I as 1, the overlap threshold of the 
polymers. The extremely large size of our chains, 
allied to a careful analysis of the size distribution 
in the dilute limit, allows us to extract an accurate 
estimate of the self-avoiding walk susceptibility ex- 
ponent 7 = 1.165 ±0.01. 

3. A satisfactory scaling collapse of the specific heat 
in both density regimes was obtained by plot- 
ting cJJj)/E 2 — 6 as 0.5 versus I. As suggested 
earliercJ, the maximum of the specific heat per 
monomer c v occurs at E — 2/8 = 4 for high con- 
centrations 4>- 

4. Chain conformations are described within numeri- 
cal accuracy by the same universal functions as for 
conventional polymers. 

5. As in the case of conventional polymers, the chains 
are swollen within the excluded volume blobs and 
Gaussian at larger distances - this has been made 
directly evident from the scaling of coil size against 
chain length within a single system. 

We believe that the present work unambiguously con- 
firms the scaling results for an idealized model of EP 
based on.±he classical behavior of conventional quenched 
polymeraLj. Accordingly, it leaves completely unan- 
swered the question of how, in some experimental sys- 
tems to which the model appears closely applicable, a 
growth exponent a ~ 1.2 is convincingly argued to 
ariseEJ. This question remains open, but clearly must 
involve physics not in the present model, which is based 
on assuming a fixed scission energy, athermal excluded 
volume interactions, and the absence of rings. 

A parallel version of the algorithm has also been im- 
plemented and used for some of the results presented 
above. It forms the basis of the natural extension of 
this work to the question of EP dynamics, which will be 
published elsewhereEa 
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Table 1. Summary ol measured static quantities for configurations with (L) > 50. Three decades of volume fractions 
between <f> — 0.001 and 4> = 0.6 have been sampled with bond energies up to E — 15. Note that Sb(4> < 0.1) = 200a and 
Lbox(4> > 0.1) = 100a. Quantities tabulated: the mean number of chains (M), the mean chain length (L), the specific heat 
per monomer c v , the mean end-to-end distance R e , the mean gyration radius R g and the average distance between chains H. 
All length scales are given in units of the lattice constant a. 
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Table 2. Variation of system size Sb for <f> — 0.5 for two high scission energies. We give the average number of chains (M), 
the mean chain length (L), the specific heat c v , the end-to-end distance R e and the radius of gyration R g . 
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